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We present preliminary results of matrix elements of four fermion operators relevant to the determination of e 
and e'/e using staggered fermions. 



1. INTRODUCTION 

To calculate the matrix elements relevant to 
CP violation in Kaon decays it is important to 
use a lattice formulation which preserves (some) 
chiral symmetry. In the case of Bk, B!f^'^\ and 

(3 /2) 

Bg , the absence of chiral symmetry leads to 
mixing with wrong chirality operators, which in 
turn leads to large discretization errors. The 
problem is far more severe for Bg^^^^ due to mix- 
ing with lower dimension operators. Two lattice 
formulations that respect at least part of the con- 
tinuum chiral transformations and hold promise 
for these calculations are domain wall/overlap 
fermions and staggered fermions. Each has its 
advantages and disadvantages. Renormalization 
of operators in the domain wall/overlap formula- 
tion is small enough that 1-loop calculations may 
be adequate, but the numerical simulations are 
~ 100 times more costly. Staggered simulations 
are very efficient, but the 1-loop renormalization 
constants for the simplest lattice transcription of 
operators are very large. The goal of this project 
is to find an improved staggered formulation for 
which perturbation theory is well-behaved. 

Here we present preliminary estimates of Bk, 
^(3/2)^ ^(3/2)^ ^(1/2) quenched lat- 

tices (16"^ x64) a.t f3 — 6.0. This numerical simula- 
tion is being done on the QCDSP supercomputer 
at Columbia University. To facilitate chiral ex- 
trapolations, we have used four values of quark 
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mass: arug = 0.01, 0.02, 0.03, 0.04. The results 
are for gauge invariant staggered operators that 
lie in a 2"* hypercube for which the 1-loop renor- 
malization constants are now known 

2. Bk 

We have used the calculation of Bk as a test 
of our programs. The estimates shown in Fig- 
ure [| agree with previous calculations and 
it is worth mentioning that the 1-loop calculations 
done independently in Ref. ||l| reproduce the re- 
sults given in Ref. The figure shows a fit us- 
ing the form suggested by chiral perturbation the- 
ory, Co + ciiaMK^ + C2{aMK)^\n{aMKy, with 
Co = 0.55(8), ci = 0.06(25), and C2 = -0.60(38). 
Physical kaons correspond to (aM^)^ ~ 0.06. 

3. i?f and i?f 

In Figures |2| and ||, we compare the tree level 
and one-loop results for B!f^^^ and B^^"^^ (de- 
fined in 1^). The 1-loop tadpole improved renor- 
malization constants for the gauge invariant op- 
erators were recently calculated in . We match 
to the continuum NDR scheme at q* ^ fi ~ 

(3/2) 

n/a. The dominant contribution to (Og ) 
comes from 2-color trace staggered operator [P x 
P][P X P]ii, which also dominates the vacuum 
saturation contribution. Thus, even though the 
renormalization constants are large, there is a 
close cancellation and the B parameter receives 
a < 10% 1-loop correction. On the other hand 
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Figure 1. Bk{h = tt/a, NDR). 



Figure 2. Bf''^\^i = TT/a,NDR). 



ipf^'^^ ) is dominated by the 1-color contraction, 
and the renormalization constants do not cancel. 
Consequently, even assuming q* — Tt/a, 1-loop 
perturbation theory is unreliable. 

These results can be compared against previ- 
ous calculations done using gauge non-invariant 
Landau gauge operators in Ref. Q|. The authors 
of Ref. 1^] found significant systematic differences 
between results obtained using smeared and un- 
smeared operators. Choosing the same value of (3 
and orMq, we find that our results lie in between. 

4. B^i^'^ 

Accurate estimates of e'/e require measure- 
ments of matrix elements of QCD and elec- 
tromagnetic penguin operators, (Og^^^^) and 
(Og ) respectively, in the K ^ nn transition. 
Since these two contribute with opposite sign, 
leading to a significant cancellation, both need 
to be measured precisely in order to test whether 
the Standard Model explains the observed size of 
Re(e7e). 

Direct calculations of the K —^ vrvr ampli- 
tudes on the lattice are difficult The sim- 
plest approach has been to assume chiral pertur- 



bation theory provides accurate relations between 
K ^ TTTT, — > TT and K ^ amplitudes 0. In 
this method, which we use, the operator Oq^^^^ 
has three types of contractions: eight, eye, and 
subtraction (to remove mixing with lower dimen- 
sion operators) [HM. Schematically, 



B, 



(1/2) 



o. 



1/2 



Vac. Sat. 



Oe^' = O^^' (Eight) 



06/'(Eye) 



In order to restrict the subtraction term to the 
single dimension four operator Ogub = {iTid — 
Ws)s75(i + {md + ms)sd, we need to work with 
degenerate s and d quarks. For — md we 
can determine the coefficient a of Osub by cal- 
culating the derivative of the K amplitude 
with respect to the strange quark mass. This 
derivative introduces two types of diagrams: (i) 
each strange quark propagator is replaced by its 
derivative, (ii) a disconnected diagram that arises 
from the differentiation of the fermion determi- 
nant. In principle, both terms can be estimated in 
the quenched theory. However, the disconnected 
contribution vanishes in the limit of rus — md- 
Hence, our calculation of a does not include the 
disconnected contributions and is thus directly 
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Figure 3. B^^^'^>{fi = ir/a, NDR). 



comparable to those of Refs. [^UgJOj, which do not 
discuss such contributions. The results are shown 
in Figure |4[ labeled as "standard" . The details 
will be presented in Ref. |p^ . 

Recently, Golterman and Pallante have pointed 
out a subtlety associated with the quenched 
approximation using (partially) quenched chiral 
perturbation theory They show that the 

usual quenched operators lead to matrix elements 
which have a different chiral expansion from that 
in the continuum. They propose an alternative 
quenched operator which does not have these 
problems. In practice this amounts to dropping 
certain contractions in the "eye" and "subtrac- 
tion" diagrams. 

We demonstrate in Figure ^ the difference in 
Bq using the standard operator and that sug- 
gested by Golterman and Pallante. The 1-loop 
renormalization constants in the NDR scheme are 
from Ref. |0] and ||l^ . Our calculations show that 
the Golterman-Pallante operator enhances by 
almost a factor of two, which in turn would sig- 
nificantly increase the value of e'/e compared to 
results given in jj^^ . Thus the quenched approx- 
imation remains the most significant drawback of 
such calculations, which we hope to address in 
the future. 
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Figure 4. BQ^^^\fj, = tt/o, NDR) using the stan- 
dard operator and Golterman-Pallante operator. 
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